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Abstract 


Rocks on the floor of Gusev crater are basalts of uniform composition and mineralogy. 
Olivine, the only mineral to have been identified or inferred from data by all instruments 
on the Spirit rover, is especially abundant in these rocks. These picritic basalts are 
similar in many respects to certain Martian meteorites (olivine -phyric shergottites). The 
olivine megacrysts in both have intermediate compositions, with modal abundances 
ranging up to 20-30%. Associated minerals in both include low-calcium and high- 
calcium pyroxenes, plagioclase of intermediate composition, iron-titanium-chromium 
oxides, and phosphate. These rocks also share minor element trends, reflected in their 
nickel-magnesium and chromium-magnesium ratios. Gusev basalts and shergottites 
appear to have formed from primitive magmas produced by melting an undepleted mantle 
at depth and erupted without significant fractionation. However, apparent differences 
between Gusev rocks and shergottites in their ages, plagioclase abundances, and volatile 
contents preclude direct correlation. Orbital determinations of global olivine distribution 
and compositions by thermal emission spectroscopy suggest that olivine-rich rocks may 
be widespread. Because weathering under acidic conditions preferentially attacks olivine 
and disguises such rocks beneath alteration rinds, picritic basalts formed from primitive 
magmas may even be a common component of the Martian crust formed during ancient 
and recent times. 

1. Introduction 

On Earth, abundant olivine in a basalt (such a rock is termed “picritic”) is usually an 
indication that the basaltic magma was mantle-derived and arrived at the surface without 



significant fractionation. Given the dominance of olivine in planetary mantles (reflecting 
its high abundance in the chondritic protoplanets from which planets were assembled), 
mantle melting must have produced olivine-bearing basalts on all the terrestrial planets. 
Olivine is invariably an early crystallizing phase in such magmas, resulting in the 
common occurrence of large olivine phenocrysts in lavas and their segregation to form 
cumulates within magma chambers. However, terrestrial experience suggests that picritic 
magmas rarely reach the surface, but instead are fractionated during transit from their 
mantle source regions. 

Basalts encountered near the Spirit Rover landing site in Gusev Crater (14.64°S 
175.36°E) contain high proportions of olivine, comparable to those in olivine-rich 
Martian basaltic meteorites (“olivine-phyric shergottites”) and the most olivine-rich Mars 
surface materials inferred from orbital spectroscopy. These rocks, analyzed by the 
Athena instruments (science instruments onboard the Mars Exploration Rovers), may 
yield new insights into the Martian mantle that partially melted to produce primitive 
magmas. Here, we expand an earlier, preliminary description of these basalts [McSween 
et a/., 2004] and explore possible petrologic interpretations. Comparison of Gusev 
picritic basalts with olivine-phyric shergottites allows a more rigorous interpretation of 
their origin, and comparison with regionally extensive, olivine-rich surface units mapped 
using thennal emission spectroscopy provides a global geologic context for olivine- 
bearing rocks in the Martian crust. 


2. Rock morphologies and distribution 



Gusev is a Noachian impact crater, 160 kilometers in diameter, located in the Aeolis 
region of Mars. The region of the floor of Gusev Crater where Spirit landed consists of a 
plains unit (previously mapped as a low-albedo smooth unit on the Gusev Crater 
Formation [Kuzmin et al., 2000] and as low-albedo materials on the Plains 
thennophysical unit [Milam et al., 2003]), thought to be a minimally graded landscape 
decorated with impact craters and dust devil tracks [Grant et al., 2004]. Rocks scattered 
about the plains grade in color from dark gray to light-toned (the latter are coated with 
variable amounts of dust) [Bell et al., 2004], Many faceted rocks have been sculpted by 
wind erosion. Some rocks are partially buried by fines, whereas others are perched on 
the surface. Most rocks are porphyritic with aphanitic groundmass and contain abundant 
vugs or vesicles, confirming their volcanic origin. The closest obvious volcanic source is 
Apollinaris Patera, a large volcano of Hesperian age located directly north of Gusev 
[Robinson and Mouginis-Mark, 1993]. 

During much of its first 100 sols, Spirit traversed the ejecta blanket of a 210 m crater 
informally named Bonneville (Figure 1). The abundance and average size of rocks 
increased with proximity to the crater. These rocks were apparently excavated from 
Bonneville, although no rock outcroppings were observed in the crater walls. Rocks are 
generally angular and fragmented, but no re-cemented breccias have been recognized. 
Three rocks (Adirondack, Humphrey, and Mazatzal) were imaged and analyzed as-is 
(i.e., the natural, dust-covered surface), after brushing with the Rock Abrasion Tool 
(RAT) to remove surface dust, and after grinding several millimeters into the rock 
interior by the RAT. Several other rocks (Mimi and Planck) encountered within the 
Bonneville ejecta blanket were analyzed only as-is. 



Following exploration of the Bonneville ejecta, Spirit crossed ejecta from an adjacent 
crater, Missoula. Then, during Sols 111-159, Spirit traversed the intercrater plains to the 
Columbia Hills. Two other rocks (Route 66 and Joshua) encountered on the intercrater 
plains were studied. Route 66 was brushed and then analyzed, and Joshua was analyzed 
as-is. As Spirit approached the Hills, several highly vesicular rocks (e.g. Dagwood, 
Figure 2e) and others resembling scoria were observed. Unfortunately, images of these 
rocks were not downlinked until after the rover had passed by, so they were not analyzed. 

The locations of all the rocks described and analyzed here are illustrated in a traverse 
map (Figure 1). True-color Panoramic Camera (Pancam) images of the rocks are shown 
in Figure 2. 

3. Petrography and physical properties of Gusev basalts 

3.1. Microscopic petrography 

Brushing and grinding using the RAT, followed by imaging using the Microscopic 
Imager (MI), allowed the discovery of alteration rinds on rocks and revealed the textures 
of relatively fresh rock interiors [ Herkenhoff et al., 2004; McSween et al., 2004]. 
Although MI images are monochromatic, overlay of high-resolution Pancam images has 
produced false-color images that reveal more petrographic details (Figure 3) [Bell et al., 
2004]. 

The ground (RATed) surfaces of Adirondack, Humphrey, and Mazatzal all exhibit 
abundant, dark gray grains (Figure 3), interpreted as olivine crystals [McSween et al., 
2004], in an aphanitic groundmass. (We will use the non-genetic term “megacrysts” to 
describe these crystals, all of which are larger than -100 pm - the resolution limit of MI 



images.) Microscopic images have sufficient resolution to permit point-counting, the 
first time such quantitative petrography has been attempted where the rocks and the 
operator are separated by millions of kilometers. Point counts of RATed rocks \McSween 
et al., 2004] indicated that megacrysts comprise approximately 9%, 25%, and 9% by 
volume in Adirondack, Humphrey, and Mazatzal, respectively. Adirondack and 
Mazatzal also contain 13 and 14 volume % patches of light-colored material, 
respectively. The light-colored material may be an alteration product of olivine, but no 
obvious replacement or other physical association of the two phases is apparent. 

Figure 4 shows a processed image of RATed Humphrey, in which dark grain pixels 
have been isolated by image processing, to make clearer the shapes and distribution of 
these grains. The grains interpreted as olivine are heterogeneously distributed, and 
abundant smaller grains, not so obvious in MI images, occur in the groundmass. None of 
the grains appears to be euhedral, although some subhedral outlines are apparent. Pixel 
counts of gray-scale images may be less subjective than optical point counting on a grid, 
and a mode of dark grains in Humphrey determined in this way is 20 volume %. 

However, MI images of Adirondack and Mazatzal do not have enough contrast to pennit 
accurate pixel counts of these crystals. 

Tiny, irregular veins of light-colored material crosscut all three rocks. These veins 
cannot be fractures filled with dust, because they are not red in Pancam images. A 
surface rind on Humphrey, revealed by a fortuitous oblique grind from the brushed 
surface into the interior of the stone (see Figure 2 of McSween et al, 2004 and Plate 12 of 
Squyres et al, 2004), does not contain the gray megacrysts seen in the rock’s interior. 

This observation suggests that weathering on Mars preferentially destroys olivine. 



Mineralogical and chemical differences between Gusev rock surfaces and interiors 
support this hypothesis, which is explained by the rapid dissolution of olivine under 
acidic conditions [ Hurowitz et al., 2005; Schroder et al., 2005], A more distinctive rind 
is visible on Mazatzal (Plate 13 of Squyres et al, 2004). The thin, dark rind underlying 
the dust cover is smooth and hard, and resembles a varnish. A sliver of this rind remains 
in the RATed image of Mazatzal in Figure 3. A light-colored vein crosscuts the varnish 
and thus must be younger. This dark varnish, too, is devoid of olivine. The alteration 
rinds on Gusev basalts have been described and interpreted elsewhere [Haskin et al., 
2005] and are not the focus of this paper. 

3.2. Grinding hardness 

The RAT monitored the energy expended during grinding. The grinding energy, 
divided by the volume of rock removed, is a function of several physical properties, but 
especially rock hardness. Because secondary alteration minerals are usually soft, 
grinding hardness is probably a proxy for degree of alteration. Grinding hardness for the 
three rocks analyzed, expressed in J/ninr , are approximately as follows: Adirondack - 
64, Humphrey - 94, and Mazatzal - 61 [Bartlett et al., 2005]. The hardness of Humphrey 
suggests that it is the least altered of these rocks, consistent with the MI observation that 
it contains the highest proportion of olivine. 

4. Olivine identification and accompanying mineralogy in Gusev basalts 

Olivine is the only mineral identified or inferred from data by all the Athena 
instruments on the Spirit Rover. The synergy provided by the Pancam, Mini Thermal 



Emission Spectrometer (Mini-TES), Mossbauer Spectrometer (MB), and Alpha Particle 
X-ray Spectrometer (APXS) allow documentation of the mineralogy of these rocks. 

4.1. Pancam multispectral observations 

Pancam 1 1 -wavelength (430 to 1009 nm) spectra of RAT holes ground into 
Adirondack, Humphrey, and Mazatzal are shown in Figure 5. These relative reflectance 
(radiance factor or “I/F”) spectra were calibrated using near-simultaneous measurements 
of the Pancam calibration target [Bell et al., 2005]. These are the lowest reflectivity, least 
red spectral units identified by Pancam at the Gusev landing site. The low reflectance 
and lack of a steep increase in reflectivity towards the red wavelengths indicate that the 
RAT holes are not contaminated with aeolian dust and imply the almost complete 
absence of ferric minerals within the shallow (—5-7 mm) interiors of these rocks. 

The reflectance maximum near 650-700 nm and the negative spectral slope from -700 
to 1000 nm are both consistent with the iron mineralogy of these rock interiors being 
dominated by ferrous silicate phases like olivine or pyroxene [e.g., Adams. 1974; Morris 
et al., 2000; Cloutis and Bell, 2003], The Pancam spectra do not reveal an unambiguous 
near-IR reflectivity minimum that could be used to provide a strong constraint on the 
specific ferrous mineralogy. However, the presence of a weak minimum between 900 
and 1000 nm in the Adirondack and Humphrey spectra could be consistent with the 
presence of low-calcium pyroxene, and may indicate a subtle difference in mineralogy 
between these two rocks and Mazatzal. In addition, the decrease in reflectance towards 
the longer wavelengths observed in all three spectra could be the high-energy wing of a 
mineral absorption band centered in the near-infrared beyond Pancam’s wavelength 
range. If so, this feature could be consistent with the presence of either olivine or high- 



calcium pyroxene, both of which exhibit absorptions in the 1000-1300 nm wavelength 
region. 

Further support for a ferrous mineral signature in these rock spectra comes from 
comparison of Pancam multispectral parameters to those derived from laboratory ferrous 
and ferric minerals by Clouds and Bell [2003] for similar multispectral measurements 
(Figure 6). Pancam data for the three RAT holes examined here plot in the same 
parameter space as olivine, high-calcium pyroxene, and the Martian meteorites. Similar 
constraints on the ferrous mineralogy can be inferred from comparison of the Pancam 
data to additional laboratory ferric and ferrous minerals studied by Morris et al. [2000]. 

While Pancam data cannot uniquely detennine the ferrous mineralogy of these rocks, 
the implied evidence for olivine and/or pyroxenes is consistent with Mini-TES and 
Mossbauer identification of these phases (see Sections 4.2 and 4.3). Furthermore, the 
lack of evidence for ferric phases, both within the average RAT hole Pancam spectra and 
in Pancam spectra from small numbers of pixels associated with vugs and veins within 
the exposed rock interiors, implies that there has been only very limited alteration of 
these rocks, even within the relatively near-surface regions exposed by the RAT. 

4.2. Mini-TES spectral observations 

Mini-TES measurements of non-dusty rocks on the plains of Gusev demonstrate that 
all the rocks exhibit similar spectral shapes (Figure 7a). The middle portions of these 
spectra (-1200-800 cm’ 1 ) are difficult to model with mineral spectra alone because of a 
high degree of atmospheric spectral contribution caused by downwelling radiance from a 
relatively warm atmosphere onto cold rocks. However, the low-wavenumber region 
(<600 cm’ 1 ) is much less affected [Ruff et al., 2005], It is in this range that a prominent 



emissivity peak centered at -425 cm' 1 is evident in all of the spectra. As shown by 
Christensen et al. [2004], this peak is best fit by olivine with a composition in the range 
of F035-6O (expressed as mol.% of the forsterite Mg 2 SiC>4 endmember). 

Figure 7b shows a Mini-TES measured emissivity spectrum and linearly deconvolved 
modeled spectral fit of a representative warm, dark disturbed soil in Gusev Crater. The 
measured spectrum is an average of 200 individual spectra and is corrected for 
downwelling atmospheric radiance and contributions of fine-grained surface dust 
[Christensen et al., 2004], The abundances of derived minerals that fit to the measured 
spectrum have been summed into mineral groups with accuracies of 5 to 10 volume %. 
The endmember spectra shown are scaled by these abundances to show how each 
contributes to the measured spectrum. The basaltic soil is dominated by approximately 
equal amounts of plagioclase (50 %) and pyroxene (20% high-calcium clinopyroxene and 
20% pigeonite) with detectable amounts of olivine (10 %). The olivine composition of 
the soils is F045 +/_ 5 10 10, which agrees with the APXS-derived normative olivine 
composition (Section 4.4). The derived bulk composition of Gusev soils also agrees with 
the compositions of basalts derived from Mars Global Surveyor (MGS) Thermal 
Emission Spectrometer (TES), which range from 45-50% plagioclase, 25-35% high- 
calcium clinopyroxene, 10-15% olivine, and 10-15% sheet silicates and/or amorphous 
phases [ Christensen et al., 2000] (Section 5.2). However, the modal mineralogy of the 
local soils does not coincide with those inferred for the rocks. 

4.3. Mossbauer observations 

Mossbauer spectra of the rocks are characterized by two Fe doublets, interpreted as 

3_l_ 

olivine and probably pyroxene, and a Fe doublet, suggested to represent nanophase 



ferric iron oxide and magnetite [Morris et al., 2004]. Ground rock interiors show the 
least Fc associated with nanophase oxide, implying that weathering is not extensive. 
Least squares fitting of the spectra allow estimates of the molar Fe /Fe(total) values for 
RATed rocks: Adirondack 0.83, Humphrey 0.84, and Mazatzal 0.90, all +0.04, similar to 
values for terrestrial basalts (0.85-0.90). 

The composition of olivine was estimated from the temperature dependence of the 
Mossbauer quadrupole splitting. The estimated composition is -Foeo [Morris et al., 
2004]; this composition assumes a linear relationship between composition and the 
splitting parameter, and no uncertainty was reported. Fo 6 o is slightly more magnesian 
than the olivine composition estimated from Mini-TES spectra or APXS nonn 
calculations (see Section 4.4). Because the Mossbauer spectrometer analyzes a deeper 
volume of material than Mini-TES or APXS, its spectra may sample a different 
proportion of olivine core compositions, which would be magnesian if the crystals are 
zoned. 

4.4. APXS measurements and norms 

Preliminary APXS analyses of Gusev rocks [Gellert et al., 2004] have recently been 
revised [Gellert et al., 2005], based on a new calibration. Sensor heads on the two rovers 
were inadvertently switched just before launch, and that belated discovery necessitated 
recalibration. Consequently, the results reported here differ slightly (about 10% relative 
for oxides of the low atomic number elements, generally within the quoted uncertainties 
of Gellert et al. [2004]) from those of the preliminary report; NaoO and MgO are lower 
by ~1 weight % absolute, and AI 2 O 3 is higher by the same amount. For other elements, 
there are no discemable differences, except those arising from re-nomalization. The 



revised APXS analyses and uncertainties for Adirondack, Humphrey, and Mazatzal 
(natural, brushed, and RATed) are tabulated by Gellert et al. [2005], 

Plots of oxides versus SO3, originally developed to correct for dust coverings on Mars 
Pathfinder rocks [ Rieder et al., 1997], are illustrated for Gusev basalts and soils in Figure 
8. Pathfinder rocks were variably coated with dust, and clean rock compositions were 
estimated by extrapolating the trends to 0.3 wt. % S (0.75 wt. % SO3), the average sulfur 
abundance (as igneous troilite) in Martian basaltic meteorites [ Bruckner et al., 2003]. 
Surface (as-is), brushed, and RATed rock compositions, which have progressively less 
SO3, are shown by small symbols in Figure 8. For most elements, the brushed 
compositions do not lie on straight lines connecting the RATed compositions to surface 
or soil compositions. Compositions of the brushed areas represent thin alteration rinds 
seen in MI imagery and are distinct from the soil compositions. However, natural rock 
surface compositions commonly plot along mixing lines between brushed and soil 
compositions. 

Analyses for Route 66, Mimi, Planck, and Joshua are also given by Gellert et al. 
[2005] and plotted in Figure 8. Soil analyses form a cluster at ~6 wt. % SO3, with a trend 
to lower SO3 (presumably representing admixture of local rock fragments) and another to 
higher SO3 (apparently representing addition of MgSCfi) seen in subsurface soils exposed 
in trenches [ Haskin et al., 2005], 

The larger symbols in Figure 8 represent end-member rock compositions determined 
by extrapolation from the brushed through the RATed values to 0.75 weight % SO3, 
equivalent to 0.83 weight % FeS. Visual inspection showed that grinding commonly did 
not remove all vestiges of the rinds, so this extrapolation corrects for any remaining rind 



in RAT holes. The extrapolated rock compositions are given in Table 1. Rock 
compositions based on the new calibration are only slightly different from those reported 
by McSween et al. [2004]; MgO and Na 2 0 are lower and AI 2 O 3 is higher. These rock 
compositions plot on the boundary separating the basalt and picrobasalt fields in the 
alkalis-silica diagram used to classify volcanic rocks (Figure 9), hence our use of the term 
“picritic basalt” to describe these rocks. 

4.5. Summary constraints on rock mineralogy 

Calculated CIPW norms for the new end-member rock compositions are given in 
Table 1. These calculations assume Fe 2+ /Fe(total) ratios for RATed rocks as detennined 
by Mossbauer measurements [Morris et al., 2004], as described above. The normative 
mineral proportions for the three basalts are illustrated in Figure 10. 

These data compare favorably with the mineral assemblage identified by Mossbauer 
spectroscopy, but the inferred relative proportions of minerals differ. Proportions of iron- 
bearing minerals (relative peak areas % [Morris et al, 2004], assumed to be volume %, 
and renormalized for non-iron-bearing minerals - plagioclase, apatite - in calculated 
nonns) for Adirondack/Humphrey/Mazatzal, respectively, are: olivine 35/35/40, 
pyroxenes 22/22/22, and oxides 15/14/9. However, it is not possible to compare directly 
the normative mineral abundances (weight %,) with Mossbauer abundances (volume %). 

As noted earlier, Mini-TES spectra for rocks cannot be deconvolved because of down- 
welling radiance effects, but olivine is confirmed. Pancam spectra are consistent with the 
presence of olivine and/or pyroxenes. Quantitative MI image analysis of dark crystals 
inferred to be olivine in Humphrey is similar to the normative olivine content, but olivine 
point counts in the other rocks are considerably lower. 



Of the minerals present, olivine has the most tightly constrained composition. 
Normative olivine compositions (illustrated by arrows) are compared with olivine 
compositions determined by Mossbauer and Mini-TES spectra for the same rocks in the 
lower part of Figure 11a. The average normative olivine composition for all three rocks 
is ~Fo 49 . Olivine abundances are compared in the lower part of Figure lib. 

5. Comparisons with martian meteorites and orbital spectroscopic data 

5.1. Olivine-phyric shergottites 

Olivine-phyric shergottites [ Goodrich , 2002] are extrusive pyroxene -plagioclase 
basalts with porphyritic textures of large olivine crystals (megacrysts). Related plutonic 
ultramafic rocks that also contain olivine are called lherzolitic shergottites, and related 
basalts without olivine megacrysts are called basaltic shergottites. Olivine-phyric 
shergottites provide the most direct analogy to the Gusev picritic basalts. 

The compositions and proportions of olivine in olivine-phyric shergottites and 
descriptive references for these meteorites are summarized in Table 2. A backscattered 
electron (BSE) image of the SAU 005 meteorite (Figure 12) shows a similar texture to 
those of Gusev basalts. The olivine megacrysts in these meteorites are strongly zoned 
with magnesian cores and ferroan rims, and their measured compositional ranges (bars in 
Figure 11a) generally overlap the average composition of olivine in Gusev basalts. This 
comparison is hampered by the fact that we do not know whether Gusev olivines are also 
zoned. If they are, as seems likely, their compositional ranges may be approximated by 
the compositions of olivines formed between the liquidus and solidus in a MELTS 
crystallization calculation (see Section 6.1). This range (F 081 - 55 ) is also illustrated by a 



gray bar labeled MELTS in Figure 1 la. To allow a more straightforward comparison, the 
calculated normative olivine compositions for olivine-phyric shergottites (calculated 
assuming bulk rock Fe /Fe(total) = 0.84, as measured for Humphrey [Morris et al, 
2004]) are illustrated by shaded circles in Figure 11a. The shergottite olivines appear to 
be more magnesian than those in Gusev basalts, but the differences are not pronounced. 
The compositions of olivines in lherzolitic shergottites (shaded box in Figure 11a) are 
also more magnesian. The volume proportions of olivine megacrysts in Gusev basalts 
and olivine-phyric shergottites also overlap (Figure lib and Table 2). These similarities 
support the hypothesis that olivine-phyric shergottites and Gusev basalts may have had a 
similar petrogenesis. What insights can be gained from this comparison? 

It is not clear whether the gray megacrysts (presumed to be olivine) in Gusev basalts 
are phenocrysts that grew from the magma or are foreign crystals added to the magma 
[McSween et al., 2004]. Likewise, there has been considerable controversy about the 
origin of the megacrysts in olivine-phyric shergottites. Some studies have focused on the 
ragged and sometimes embayed grain boundaries, the presence of occasional 
polymineralic grains, and the compositions of olivine cores that are generally too iron- 
rich to have been in equilibrium with the bulk-rock compositions. Such studies have also 
noted that ferroan rims occur on megacrysts where they are in contact with groundmass, 
suggesting disequilibrium with the enclosing magma. Based on such observations, the 
olivines have been suggested to be xenocrysts (foreign crystals incorporated into the 
magma) [Steele and Smith, 1982; McSween and Jarosewich, 1983, Wadhwa et al., 2001, 
Koizumi et al, 2004] or cumulates (crystals accumulated and concentrated from a large 
batch of magma) [Barr at et al., 2002, Koizumi et al., 2004], In these models, olivine has 



also been suggested to have been partly resorbed and assimilated by the host magmas. 
The high thennal energy required to assimilate olivine crystals poses a considerable 
problem, and several alternatives have been proposed. Wadhwa et al. [1994] suggested 
that the megacrysts might have been introduced by mixing an olivine phenocryst-bearing 
magma with a basaltic shergottite magma. Mittlefehldt et al. [1999] favored the idea that 
the megacrysts were residual crystals from impact melting, and Folco et al. [2000] 
proposed that they were restites from partial melting of lherzolitic shergottites. 

Recent recovery of new olivine-phyric shergottites with distinct textures and more 
magnesian olivine core compositions supports another interpretation. Some of these 
(Dhofar 019 and Yamato 980459) contain euhdral olivine megacrysts having magnesian 
core compositions (F 083 - 86 ) that are in equilibrium with the bulk-rock compositions 
[Taylor et al., 2002; Greshake et al., 2004; Mikouchi et al., 2004], Crystal size 
distribution (CSD) analyses of olivines for the most part produce linear arrays [ Taylor et 
al., 2002; Goodrich, 2003; Greshake et al., 2004], supporting the hypothesis that most 
olivines formed by continuous cooling without interruption. These data support the 
hypothesis that the olivine megacrysts are phenocrysts. 

A complicating factor is that even Dhofar 019 and Yamato 980459 contain a few 
corroded, Fe-rich olivines, and their CSD patterns deviate from the linear trend at the 
largest grain sizes [ Goodrich , 2003; Greshake et al., 2004]. These unusual grains are 
possibly xenocrysts or cumulates, but they comprise only a small portion of the olivine 
megacryst population in each meteorite. A further complication is that, for most olivine- 
phyric shergottites, the earliest crystallizing olivine phenocrysts (those with compositions 
in equilibrium with the bulk-rock magma composition) are missing, apparently removed 



by crystal fractionation [ Goodrich , 2003]. Taking all the data into account, the most 
plausible model is that the bulk of olivine crystals in olivine -phyric shergottites are 
phenocrysts. In many of these meteorites, the earliest-formed phenocrysts have been lost 
and some small proportion of xenocrysts or cumulates have been added. However, the 
most magnesian olivine -phyric shergottites like Yamato 980459 appear to represent 
liquid compositions, or nearly so. 

As expected, olivine-phyric shergottites have higher contents of magnesium and 
nickel than do basaltic shergottites (Figure 13a). Given the model above, this correlation 
can be interpreted as reflecting the higher magnesium and nickel abundances in primitive 
magmas, with the trend representing an olivine-control fractionation line. The Gusev 
basalts Humphrey, Adirondack, and Mazatzal plot near the olivine-phyric shergottites in 
Figure 13a. All the Martian rocks have higher magnesium contents at a given nickel 
content than terrestrial mafic and ultramafic rocks (Figure 13a). Gusev basalts also have 
magnesium-chromium ratios that are similar to olivine-phyric shergottites (Figure 13b). 

In this case, chromium occurs dominantly in associated spinel (also an early-crystallizing 
phase, see Section 6.1) rather than being hosted in olivine. 

The tight clustering of Gusev basalt compositions may provide a further argument that 
the olivines they contain are phenocrysts. Incorporation of xenocrysts or cumulates 
would likely be a random process, producing variable proportions of megacrysts and 
varying bulk rock chemistry. Despite apparent differences in modal olivine contents in 
Adirondack, Humphrey, and Mazatzal, their chemical compositions are nearly uniform. 
5.2. Orbital spectroscopy and regional/global context 



The Thermal Emission Spectrometer (TES) on Mars Global Surveyor has provided 
unprecedented insights into the composition of broad areas of the Martian surface. A 
detailed study of atmospherically corrected spectra [Smith et al., 2000] in Cimmeria Terra 
by Christensen et al. [2000] identified basaltic surface compositions dominated by 
plagioclase (45% and 53%) and high-calcium pyroxene (26% and 19%) with detectable 
amounts of olivine (12%). Further analyses [ Hoefen et al., 2003; Hamilton and 
Christensen, 2005] of TES spectra using several intermediate olivine compositions 
suggested that olivine is widespread and abundant in Nili Fossae, a volcanic terrain 
northeast of Syrtis Major. Mapped olivine compositions in Nili Fossae overlap olivine 
compositions in Gusev basalts, but olivines in Syrtis Major are more iron-rich (Figure 
1 la). The olivine composition becomes progressively more iron-rich from west to east 
across the 30,000 km Nili Fossae area, and linear deconvolutions indicate an olivine 
abundance of -30%. A study by Hamilton et al. [2003] also found high abundances of 
olivine in Nili Fossae and local concentrations including but not limited to Gangis 
Chasma, Aurorae Planum, the Argyre and Hellas basin rims, and Eos Chasma. Rogers et 
al. [2005] also found several olivine-bearing bedrock units, containing 25% olivine of 
composition Fo6o, in Ares Vallis using data from TES and the Thermal Emission Imaging 
System (THEMIS). These authors noted that olivine -phyric shergottites provide the 
closest compositional match for the Ares Vallis rocks. Their exposure patterns suggest 
formation as several distinct flows. 

Figure 14a shows a THEMIS daytime infrared mosaic of Gusev Crater. 

Superimposed on the mosaic are the MER Spirit landing ellipse (-83 km long by -10 km 
wide) and a TES orbital track (ock 5367, ick 1600-1603) with different colors showing 



derived surface temperatures (3x6 km spatial resolution). Differential heating of slopes 
produces a clear picture of the crater rim, Ma’adim Vallis, and mesas in the south. The 
floor of Gusev Crater, however, displays temperature variations that are independent of 
topography. Two prominent streaks trending SSE are warmer in the day than the rest of 
the crater floor and correspond to relatively low albedo, dust-free streaks that are evident 
in visible images. Their warmer temperatures are due to the differential heating that 
result from albedo contrast. 

An average atmospherically-corrected emissivity spectrum of six of the wannest TES 
pixels (294-296K) outlined in the THEMIS mosaic is shown in Figure 14b with a linearly 
deconvolved modeled spectral fit. The abundances of derived minerals that fit the 
measured spectrum have been summed into mineral groups with accuracies of 5-10 
volume %. The endmember spectra shown are scaled by these abundances to show how 
each contributes to the measured spectrum. The composition of the warm, low-albedo 
streak is dominated by plagioclase (45 %) and clinopyroxene (30 %) with 15% olivine. 
This bulk mineralogy agrees well with the Mini-TES derived composition (Section 4.2) 
of dark disturbed soils (50 % plagioclase, 40 % clinopyroxene, and 10 % olivine) and 
indicates a regional similarity in basaltic compositions in Gusev Crater. 

Figure 15 shows a new global map of olivine abundances on Mars derived from 
atmospherically-corrected TES data, overlain on a color mosaic. Recent work [ Hamilton 
and Schneider, 2005] has led to four additional intennediate to ferroan olivine 
composition (Fo68, Fo6o, F035, and F010) spectral endmembers that can be used to detect 
and model olivine abundances more accurately. TES emissivity spectra were binned and 
averaged at 4 pixels/degree (15 km/pixel) [ Bandfield et al., 2002] and linearly 



deconvolved from 1301-825 cm' 1 and 508-233 cm' 1 using a spectral endmember set that 
includes a broad range of igneous and sedimentary minerals, Martian atmospheric dust 
and water ice [Bandfield et al., 2000b], and epf (emission phase function)-derived 
Martian surface dust [Bandfield and Smith, 2003]. The new TES olivine map shows 
detectable amounts of olivine (10-15 volume %) in broad near-equatorial regions and 
corresponds well with the highest mapped abundances of the global TES basalt unit 
[Bandfield et al., 2000a], This result is consistent with previous work by Christensen et 
al. [2000] which modeled 12% olivine for basaltic surface compositions. The global 
olivine map also resolves higher olivine abundances (35-50 volume %) detected in the 
Nili Fossae region [Hoefen et al., 2003; Hamilton and Christensen, 2005] but does not 
resolve olivine abundance in Gusev Crater because of the decreased spatial resolution (15 
km/pixel) of the binned dataset and high amount of surface dust. 

Results from individual TES spectra of Gusev Crater and globally mapped olivine 
abundances indicate that large expanses of the Martian surface are characterized by 
olivine-bearing basalts and/or soils with similar compositions. 

6. Discussion 

6.1. Do Gusev picritic basalts represent primitive magmas? 

We have already seen that it is difficult to detennine whether the olivines in Gusev 
basalts are phenocrysts or were added to magmas as xenocrysts or cumulate material. 

The 1-atm crystallization sequence of a magma having the bulk composition of 
Humphrey (chosen because it appears to be least altered, based on MI observations and 
grinding hardness) has been estimated using the MELTS program [Ghiorso and Sack, 
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1995]. The measured Fe /Fe(total) value for Humphrey corresponds to an oxidation 
state close to the quartz-fayalite-magnetite buffer (QFM), so we assumed that oxidation 
state for this calculation. The results (Figure 16) indicate early crystallization of 
chromium spinel and then olivine, followed by the near-simultaneous appearance of 
orthopyroxene, clinopyroxene, plagioclase, and whitlockite. One possible interpretation 
of the early appearance of olivine and spinel is that both phases accumulated in the 
magma or were incorporated as xenocrysts. 

However, high-pressure experiments reported by Monders et al. [2005] indicate that 
the Gusev picritic basalt magma composition is multiply saturated with olivine, spinel, 
and orthopyroxene at ~1 1 GPa, corresponding to a depth in the Martian mantle of -75 
km. This experiment offers strong support for the idea that Gusev basalts represent a 
primitive magma composition, formed by -20% melting of an undepleted mantle source. 
Other subsolidus experiments are consistent with a Martian mantle of olivine + 
orthopyroxene + spinel at this depth [Bertka and Holloway , 1994], This magma may 
have experienced little or no fractionation during ascent and emplacement, but instead 
underwent equilibrium crystallization at or near the surface. 

6.2. Could Gusev basalts be olivine-phyric shergottites? 

We previously noted similarities in the modal abundance and composition of olivines 
in Gusev basalts and olivine-phyric shergottites, as well as in their bulk-rock 
nickel/magnesium and chromium/magnesium ratios. There is yet another intriguing 
similarity that should be noted. 

Preliminary high-pressure melting experiments on Yamato 980459, the most 
magnesian olivine-phyric shergottite [ Greshake et al, 2004; Mikouchi et al., 2004], 



indicate co-saturation with olivine and pyroxene at 1 1-12 GPa [Dalton et al., 2005]. 
Olivine compositions in the experiment (Fo86) match those of the natural olivine cores. 
Thus, Gusev basalts and olivine-phyric shergottites appear to share origins as primitive 
magmas formed at comparable mantle depths. 

Nonetheless, there are some significant compositional differences between these 
rocks. Shergottites have lower AI2O3 and Na?0 (reflected in lower plagioclase 
abundance) and higher Si02 (Figure 9) than Gusev basalts. On a Ca/Si versus Mg/Si 
diagram used for classification of Martian meteorites [ Oura et al., 2003], Gusev basalts 
do not plot within any of the olivine-phyric shergottite field, nor any other shergottite 
field (Figure 17). Also, the Gusev basalts commonly contain vesicles and vugs (Figure 
2e), testifying to much higher contents of volatiles on eruption, relative to the olivine- 
phyric shergottites which do not contain vesicles and were essentially dry. 

More importantly, there are apparent age differences. The measured radiometric ages 
of olivine-phyric shergottites range from 175 to 575 Ma [Nyquist et al. 200 1 ; Borg et al. 
2001], corresponding to the Late Amazonian period. Although no specific sources for 
these meteorites have been determined, they probably were derived from impact craters 
in the Tharsis or Elyssium volcanic provinces, based on the relatively young crater 
counting ages for these areas [McSween, 2002] and their lack of detection in spectral 
studies of martian dark regions [Hamilton et al., 2003]. Martinez-Alonso et al. [2005] 
suggested that Gusev basalts likely erupted from Apollinaris Patera, based on the 
presence of flows on the volcano’s southern fla nk and the existence of a favorable path 
through the northwest rim of Gusev (now blocked by a younger crater). Such lavas 
would have to travel -250 km up a significant topographic gradient, so this source is 



problematic. However, the calculated viscosity for Gusev basalt compositions [ Greeley 
et al., 2005] suggests these flows may have been more fluid than terrestrial lavas and 
comparable to mare lavas on the Moon. There are no obvious volcanic edifices nearer to 
Gusev, so we cannot specify the source of these flows. Whatever their source, the 
volcanic units surrounding Apollinaris and the plains unit within Gusev have Hesperian 
[Milam et al., 2003] or Early Amazonian to Hesperian model ages [Kuzmin et al., 2000], 
signficantly older than olivine-phyric shergottites. Morphological data and crater counts 
from orbital imagery suggest that lavas flooded Gusev at -3.65 Ga, postdating other 
materials such as putative sediments emplaced by Ma’adim Vallis [Greeley et al., 2005]. 

The differences in composition and age appear to rule out any direct petrogenetic 
relationship between Gusev picritic basalts and shergottites. However, it seems clear that 
undepleted martian mantle sources have produced primitive magmas over much of post- 
Noachian time. 

6.3. Are rocks formed from primitive magmas widespread on Mars? 

There are conflicting views about whether the Martian crust is dominated by basalt 
[Wyatt and McSween, 2002; McSween et al., 2003] or contains significant amounts of 
andesitic rocks [Bandfield et al., 2000a; Hamilton et al., 2001], Even the basalt-only 
advocates, however, have not argued that primitive basalts might be abundant. Primitive, 
mantle-derived magmas are uncommon on Earth, because most magmas experience 
fractionation and sometimes contamination with crust en route to the surface. In 
identifying primitive magmas, even olivine -rich basaltic rocks can be misleading, 
because they may be fractionated basalts containing cumulus olivine. Yet, the very 
limited surface exploration of Mars by rovers and the probable biased sampling of the 



surface by Martian meteorites may have both found instances of rocks formed by 
solidification of primitive, olivine-rich magmas. The identification of other olivine- 
bearing volcanic materials from orbital thennal emission spectra further suggests that 
primitive magmas may not be unusual. Picritic basalts may even be common, because 
Martian weathering under acidic conditions preferentially attacks olivine [ Hurowitz et al., 
2005; Schroder et al., 2005] and disguises such rocks under alteration rinds, as observed 
at Gusev. 

7. Summary 

Volcanic rocks excavated from craters on the floor of Gusev Crater are unifonn in 
composition. These picritic basalts share many similarities with olivine -phyric 
shergottites, including: 

• Intermediate to ferroan olivine compositions - Gusev basalts contain olivines with 
average compositions in the range F 060 - 40 ; shergottites contain normally zoned 
olivines ranging from Fo 84 to F025 (zoning within individual meteorites is more 
restricted) with normative olivine compositions of ~Fo 65 . 

• Olivine modal abundances ranging up to 20-30 volume %. 

• Coexisting minerals that include low-calcium and high-calcium pyroxenes, 
plagioclase of intennediate composition, iron-titanium-chromium oxides, and 
phosphate. 

• Porphyritic textures produced by olivine megacrysts. 

• Similar nickel-magnesium and chromium-magnesium systematics. 



• Major element bulk-rock compositions that are multiply saturated with olivine, 
orthopyroxene, and spinel at mantle pressures, suggesting they are primitive 
magmas. 

A number of differences are also apparent: 

• Textures - Gusev basalts have abundant vesicles and vugs, suggesting eruption 
with high volatile contents. 

• Chemical compositions - Gusev basalts have higher AI 2 O 3 and Na 2 0 contents, 
reflected in higher plagioclase abundances than in shergottites, and they do not 
plot with shergottites on a Ca/Si versus Mg/Si classification diagram. 

• Age - Gusev basalts appear to be Hesperian or Early Amazonian, whereas olivine- 
phyric shergottites have radiometric ages that are significantly younger (Late 
Amazonian). 

Examination of orbital thermal emission spectra provides a global 
context for olivine-rich basaltic rocks. Such rocks, containing intennediate to ferroan 
olivine with abundances up to 30%, appear to be widely distributed within the ancient 
highlands of Mars. The existence of younger picritic basalts, although clearly 
demonstrated by olivine -phyric shergottites, has not been confirmed by spectroscopy, 
probably because young volcanic terrains in Tharsis and Elyssium (likely sources for 
shergottites) are obscured by dust cover. Nevertheless, primitive olivine-rich basalts of 
various ages have now been suggested by rover exploration at Gusev Crater, studies of 
Martian meteorites, and orbital remote sensing. Given that weathering processes under 
acidic conditions preferentially destroy olivine and obscure primary rock mineralogy 



under alteration rinds, olivine-rich basalts formed from primitive magmas may constitute 
a significant portion of the Martian crust. 
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Figure Captions 


Figure 1. Traverse map showing locations of analyzed rocks described in this paper. The 
numbers identify the sols during which each RATed rock was analyzed by APXS. 

Figure 2. True-color Pancam mosaic images [(e) is approximate color rendering because 
some filters are missing] of basaltic rocks in Gusev Crater. All rocks are named. The 
scale bar in each image is 10 cm. 

Figure 3. MI images (colorized by Pancam overlays) of RAT holes in rocks - from top to 
bottom, Adirondack, Mazatzal, Humphrey. Humphrey is illuminated by direct sunlight 
from the top, whereas images of the other rocks were acquired while they were in 
shadow. 

Figure 4. (a) Dark crystals in Humphrey, and (b) a processed image illustrating the 
distribution and shapes of these grains (black). 

Figure 5. Pancam spectra of RAT holes ground into Adirondack (Sol 35, imaging 
sequence P2578), Humphrey (Sol 60, imaging sequence P2597), and Mazatzal (Sol 82, 
imaging sequence P2590). Errors bars represent the standard deviations of several 
hundred pixels within selected regions-of-interest. 

Figure 6. Pancam RAT hole multispectral parameters compared to laboratory ferric and 
ferrous minerals characterized by Cloutis and Bell [2003], The laboratory measurements 
were resampled to the Imager for Mars Pathfinder filter bandpasses, which are essentially 
the same as the Pancam bandpasses for the parameters plotted here. The black rectangles 
connected by the line are parameters derived from ferrous-ferric mixtures. 

Figure 7. Mini-TES spectra of Gusev rocks and soils, (a) Spectra of non-dusty rocks are 
similar, demonstrating that all the rocks on the plains have similar compositions. The 



absorption band at ~500 cm" 1 results from olivine, (b) Deconvolved spectra of average 
disturbed soils on the Gusev plains probably reflect the mineralogy of the rocks. CO 2 
marks the position of the atmospheric CO? band. 

Figure 8. APXS analyses of rocks and soils [ Gellert et al, 2005], oxides plotted versus 
SO3. Clean rock compositions (large symbols) are derived by extrapolating from brushed 
compositions through RATed compositions to 0.3 wt % sulfur (0.75 wt % SO3). Route 
66 (RATed composition) is distinguished from Mimi, Planck, and Joshua (unbrushed 
compositions). 

Figure 9. Alkalis versus silica diagram used for volcanic rock classification. Gusev 
basalts plot on the boundary between basalt and picrobasalt, distinct from Martian 
basaltic meteorites (shergottites and nakhlites, references in Table 2) and Mars Pathfinder 
rocks [ Bruckner et al., 2003; Foley et al, 2003]. The cloud of points represents global 
surface compositions estimated from TES spectra, as described by McSween et al. 

[2003], 

Figure 10. Calculated normative mineralogy for Adirondack, Humphrey, and Mazatzal 
(Table 1). Mineral abundances (wt %) are given, in order, for each of the analyzed rocks. 
Figure 11. (a) Comparison of olivine compositions (mol % forsterite MgiSiCfi) for 
Adirondack, Humphrey, and Mazatzal measured by Mini-TES and Mossbauer (MB) 
spectroscopy and APXS normative calculations with electron probe analyses for olivines 
in olivine-phyric and lherzolitic shergottites (references in Table 2). The shaded circles 
on olivine-phyric shergottite range bars are the normative olivine compositions for those 
meteorites (Table 2). The bar labeled MELTS represents the range of olivine 
compositions predicted to crystallize from Gusev basalt magmas (see text). Also shown 



are the ranges of olivine compositions determined from orbital TES spectra [Hoefen et 
al., 2003; Hamilton et al., 2003], (b) Comparison of olivine abundances in Adirondack, 
Humphrey, and Mazatzal with point-counts of olivine -phyric shergottites (Table 2; 
ranges reflect counts of multiple thin sections) and olivine abundances estimated from 
TES spectra [Hoefen et al, 2003]. 

Figure 12. Backscattered electron (BSE) image of olivine-phyric shergottite SAU 005 
showing olivine megacrysts in a groundmass of pyroxenes and plagioclase. Figure is 5.5 
mm across [after Goodrich, 2003], 

Figure 13. (a) Nickel versus magnesium in Gusev basalts, compared to bulk shergottites 
and terrestrial basalts and ultramafic rocks [modified from Zipfel et al., 2000]. (b) 
Chromium versus magnesium in Gusev basalts, compared to bulk shergottites [modified 
from Weinke and Dreibus 1988]. 

Figure 14. (a) THEMIS daytime infrared mosaic of Gusev Crater with superimposed 
MER Spirit landing ellipse and TES orbital track showing derived surface temperatures. 
Temperature variations on the floor of Gusev are due to differential heating that result 
from albedo contrast, (b) Average atmospherically-corrected emissivity spectrum of six 
of the wannest TES pixels outlined in the THEMIS mosaic with a linearly deconvolved 
modeled spectral fit. The composition of the warm, low-albedo streak in Gusev Crater 
agrees well with Mini-TES derived compositions of dark disturbed soils. 

Figure 15. Olivine abundance map of Mars, showing the widespread global distribution 
of this mineral. 

Figure 16. Calculated equilibrium crystallization sequence for a magma having the bulk 
composition of Humphrey at an oxidation state corresponding to QFM, using the MELTS 



program. The compositions of olivine (mol % forsterite) and plagioclase (mol % 
anorthite) during magma crystallization are given on the respective bars, and the total 
amounts of each phase crystallized (g of solid produced per 100 g of magma) are given at 
the bottom of the figure. 

Figure 17. Chemical classification diagram for Martian meteorites [after Oura et al, 
2003], illustrating a compositional difference between Gusev basalts and olivine -phyric 
shergottites. 
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Table 1. Extrapolated end-member chemical compositions 
and calculated norms for Gusev basalts. 



Adirondack 

Humphrey 

Mazatzal 

Oxides (wt %>) 



Si0 2 

45.30 

45.85 

45.55 

Ti0 2 

0.49 

0.58 

0.57 

A1 2 0 3 

10.42 

10.40 

10.72 

Fe 2 0 3 

3.63 

3.36 

2.11 

Cr 2 0 3 

0.62 

0.67 

0.59 

FeO 

15.71 

15.67 

16.83 

MnO 

0.42 

0.42 

0.43 

MgO 

11.90 

10.67 

10.34 

CaO 

7.76 

8.15 

8.23 

Na 2 0 

2.09 

2.35 

2.62 

k 2 o 

0.03 

0.09 

0.11 

P2O5 

0.54 

0.59 

0.63 

FeS 

0.83 

0.83 

0.83 

Total 

99.81 

99.79 

99.63 

Norms (wt %) 



plug 

36.8 (An 52 ) 

38.0 (An 46 ) 

40.0 (An 43 ) 

or 

0.18 

0.53 

0.65 

ab 

17.68 

19.88 

22.17 

an 

18.97 

17.57 

17.17 

di 

13.17 

15.70 

16.23 

by 

18.86 

16.60 

7.42 

ol 

21.7 (Fo 52 ) 

20.2 (F 049 ) 

28.6 (Fo 45 ) 

fa 

11.18 

9.86 

12.89 

fa 

10.51 

10.33 

15.73 

mt 

5.26 

4.87 

3.06 

cm 

0.91 

0.99 

0.87 

il 

0.93 

1.10 

1.08 

ap 

1.28 

1.40 

1.49 


Rock end-member compositions were calculated by subtracting from the RAT-abraded 
APXS compositions enough of the brushed APXS analyses (31%, 21%, 13% for 
Adirondack, Humphrey, and Mazatzal, respectively) to reduce residual S to 0.3 wt %. 
After conversion of S to FeS (0.83 wt %), the remaining iron was partitioned between 
FeO and Fe 2 0 3 on the basis of Mossbauer measurements of molar Fe /Fe(total) of 0.17, 
0.16, and 0.10 for the three rocks respectively [Morris et al, 2004]. Compositions were 
recast into CIPW normative mineralogy after removing FeS. Norm abbreviations are: 



plag plagioclase, or orthoclase, ab albite, an anorthite, di diopside, hy hypersthene, ol 
olivine, fo forstcritc, /h fayalite, mt magnetite, cm chromite, il ilmenite, ap apatite. 



Table 2. Properties of olivine in olivine-phyric shergotttites. 


Meteorite 

Compositional 
Range (% Fo) 

Abundance 
(vol %) 

Normative* 
composition (% Fo) 

References 

Dho 019 

72-25 

7-12 

57 

(1) 

EETA 79001 

81-52 

7-13 

62 

(2), (3) 

DaG 476+ 

79-62 

10-24 

66 

(4-7) 

DaG 489+ 

79-59 

18-20 

66 

(6), (8) 

SAU 005 

74-62 

21-29 

65 

(9) 

NWA 1068 

72-42 

21 

# 

(10) 

Y-980459 

84-31 

9-16 

65 

(11-14) 


*Norm calculations assumesFe" /Fe(total) = 0.84, as in Humphrey [ Morris et al. 2004]. 
+Probably paired, along with DaG 670, 735, 876, 975. 

#No silica analysis reported, so norm cannot be calculated. 

References: (1) Taylor et al. [2002], (2) Steele and Smith [1982], (3) McSween and 
Jarosewich [1983], (4) Zipfel et al. [2000], (5) Mikouchi et al. [2001], (6) Wadhwa et al. 
[2001], (7) Koizumi et al. [2004], (8) Folco et al. [2000], (9) Goodrich [2003], (10) 

Barr at et al. [2002], (11) Gres hake et al. [2004], (12) Ikeda [2004], (13) Mikouchi et al. 
[2004], (14) Shirai and Ebihara [2004]. 





